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(54) Assisted ventilation to match patient respiratory need 



(57) A method for calculating the instantaneous 
inspired volume of a subject during ventilatory assist- 
ance, comprises the steps of: 

estimating respiratory airflow, 
calculating an index equal to the fuzzy extent to 
which any leak is suddenly changing, and 
calculating instantaneous inspired volume as a 
function of both the estimated respiratory airflow 
and said index. 
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Description 

Field of the Invention 

5 [0001] The invention relates to methods and apparatus for the provision of ventilatory assistance matched to a sub- 
ject's respiratory need. The ventilatory assistance can be for a subject who is either spontaneously or non-spontane- 
ously breathing, or moves between these breathing states. The invention is especially suitable for, but not limited to, 
spontaneously breathing human subjects requiring longterm ventilatory assistance, particularly during sleep. 

10 Background of the Invention 

[0002] Subjects with severe lung disease, chest wall disease, neuromuscular disease, or diseases of respiratory 
control may require in-hospital mechanical ventilatory assistance, followed by longterm home mechanical ventilatory 
assistance, particularly during sleep. The ventilator delivers air or air enriched with oxygen to the subject, via an inter- 

15 face such as a nosemask, at a pressure that is higher during inspiration and lower during expiration. 

[0003] In the awake state, and while waiting to go to sleep, the subject's ventilatory pattern is variable in rate and 
depth. Most known ventilatory devices do not accurately match the amplitude and phase of mask pressure to the sub- 
ject's spontaneous efforts, leading to discomfort or panic. Larger amounts of asynchrony also reduce the efficiency of 
the device. During sleep, there are changes in the neural control of breathing as well as the mechanics of the subject's 

20 airways, respiratory muscles and chest wall, leading to a need for substantially increased ventilatory support. There- 
fore, unless the device can automatically adjust the degree of support, the amplitude of delivered pressure will either 
be inadequate during sleep, or must be excessive in the awake state. This is particularly important in subjects with 
abnormalities of respiratory control, for example central hypoventilation syndromes, such as Obesity Hypoventilation 
Syndrome, where there is inadequate chemoreceptor drive, or Cheyne Stokes breathing such as in patients with severe 

25 cardiac failure or after a stroke, where there is excessive or unstable chemoreceptor drive. 

Furthermore, during sleep there are inevitably large leaks between mask and subject, or at the subjects mouth if this is 
left free. Such leaks worsen the error in matching the phase and magnitude of the machine's effort to the subject's 
needs, and, in the case of mouth leak, reduce the effectiveness of the ventilatory support. 
[0004] Ideally a ventilatory assistance device should simultaneously address the following goals: 

30 

(i) While the subject is awake and making substantial ventilatory efforts, the delivered assistance should be closely 
matched in phase with the patient's efforts. 

(ii) The machine should automatically adjust the degree of assistance to maintain at least a specified minimum ven- 
tilation, without relying on the integrity of the subject's chemoreflexes. 

35 (iii) It should continue to work correctly in the presence of large leaks. 

[0005] Most simple home ventilators either deliver a fixed volume, or cycle between two fixed pressures. They do 
so either at a fixed rate, or are triggered by the patient's spontaneous efforts, or both. All such simple devices foil to 
meet goal (ii) of adjusting the degree of assistance to maintain at least a given ventilation. They also largely fail to meet 
40 goal (i) of closely matching the subjects respiratory phase: timed devices make no attempt to synchronize with the sub- 
jects efforts; triggered devices attempt to synchronize the start and end of the breath with the subject's efforts, but 
make no attempt to tailor the instantaneous pressure during a breath to the subject's efforts. Furthermore, the triggering 
tends to fail in the presence of leaks, thus failing goal (iii). 

[0006] The broad family of servo- ventilators known for at least 20 years measure ventilation and adjust the degree 
45 of assistance to maintain ventilation at or above a specified level, thus meeting goal (ii), but they still fail to meet goal (i) 
of closely matching the phase of the subject's spontaneous efforts, for the reasons given above. No attempt is made to 
meet goal (iii). 

[0007] Proportional assistist ventilation (PAV), as taught by Dr Magdy Younes, for example in Principles and Prac- 
tice of Mechanical Ventilation, chapter 1 5, aims to tailor the pressure vs time profile within a breath to partially or com- 
50 pletely unload the subject's resistive and elastic work, while minimizing the airway pressure required to achieve the 
desired ventilation. During the inspiratory half-cycle, the administered pressure takes the form: 

P(t) = P 0 + R.f RESP (t) + E.V(t) 

55 where R is a percentage of the resistance of the airway, f R ESP(0 is tne instantaneous respiratory airflow at time t, E is 
a percentage of the elastance of lung and chest wall, and V(t) is the volume inspired since the start of inspiration to the 
present moment. During the expiratory half-cycle, V(t) is taken as zero, to produce passive expiration. 
[0008] An advantage of proportional assist ventilation during spontaneous breathing is that the degree of assist- 
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ance is automatically adjusted to suit the subject's immediate needs and their pattern of breathing, and is therefore 
comfortable in the spontaneously breathing subject. However, there are at least two important disadvantages. Firstly, 
V(t) is calculated as the integral of flow with respect to time since the start of inspiration. A disadvantage of calculating 
V(t) in this way is that, in the presence of leaks, the integral of the flow through the leak will be included in V(t), resulting 

5 in an overestimate of V(t) ( in turn resulting in a runaway increase In the administered pressure. This can be distressing 
to the subject. Secondly, PAV relies on the subject's chemoreceptor reflexes to monitor the composition of the arterial 
blood, and thereby set the level of spontaneous effort. The PAV device then amplifies this spontaneous effort In sub- 
jects with abnormal chemoreceptor reflexes, the spontaneous efforts may either cease entirely, or become unrelated to 
the composition of the arterial blood, and amplification of these efforts will yield inadequate ventilation. In patients with 

10 existing Cheyne Stokes breathing during sleep, PAV will by design amplify the subject's waxing and waning breathing 
efforts, and actually make matters worse by exaggerating the disturbance. Thus RAV substantially meets goal (i) of pro- 
viding assistance in phase with the subject's spontaneous ventilation, but cannot meet goal (ii) of adjusting the depth of 
assistance if the subject has inadequate chemoreflexes, and does not satisfactorily meet goal (iii). 
[0009] Thus there are known devices that meet each of the above goals, but there is no device that meets all the 

is goals simultaneously. Additionally, it is desirable to provide improvements over the prior art directed to any one of the 
stated goals. 

[0010] Therefore, the present invention seeks to achieve, at least partially, one or more of the following: 

(i) to match the phase and degree of assistance to the subject's spontaneous efforts when ventilation is well above 
20 a target ventilation, 

(ii) to automatically adjust the degree of assistance to maintain at least a specified minimum average ventilation 
without relying on the integrity of the subject's chemoreflexes and to damp out instabilities in the spontaneous ver- 
itilatory efforts, such as Cheyne Stokes breathing. 

(iii) to provide some immunity to the effects of sudden leaks. 

25 

Disclosure of the Invention 

[0011] In what follows, a fuzzy membership function is taken as returning a value between zero and unity, fuzzy 
intersection A AND B is the smaller of A and B, fuzzy union A OR B is the larger of A and B t and fuzzy negation NOT 
30 A is 1 - A 

[0012] The invention discloses the determination of the instantaneous phase in the respiratory cycle as a continu- 
ous variable. 

[0013] The invention further discloses a method for calculating the instantaneous phase in the respiratory cycle 
including at least the steps of determining that if the instantaneous airflow is small and increasing fast, then it is dose 
35 to start of inspiration, if the instantaneous airflow is large and steady, then it is dose to mid-inspiration, if the instanta- 
neous airflow is small and decreasing fast, then it is close to mid-expiration, if the instantaneous airflow is zero and 
steady, then it is during an end-expiratory pause, and airflow conditions intermediate between the above are assodated 
with correspondingly intermediate phases. 

[0014] The invention further discloses a method for determining the instantaneous phase in the respiratory cycle 
40 as a continuous variable from 0 to 1 revolution, the method comprising the steps of: 

selecting at least two identifiable features F N of a prototype flow-vs-time waveform f(t) similar to an expected res- 
piratory flow-vs-time waveform, and for each said feature: 

determining by inspection the phase $ N in the respiratory cycle for said feature, assigning a weight W N to said 
45 phase, 

defining a "magnitude" fuzzy set M N whose membership function is a function of respiratory airflow, and a "rate of 
change" fuzzy set C N , whose membership function is a function of the time derivative of respiratory airflow, chosen 
such that the fuzzy intersection M N and C n will be larger for points on the generalized prototype respiratory wave- 
form whose phase is closer to the said feature F N than for points closer to all other selected features, 
so setting the fuzzy inference rule R N for the selected feature F^ to be: If flow is M N and rate of change of flow is C N 

then phase = $ N , with weight W N . 
measuring leak-corrected respiratory airflow, 

for each feature F N calculating fuzzy membership in fuzzy sets M N and C Nt 

for each feature F N applying fuzzy inference rule to determine the fuzzy extent 

55 

Y N = M N AND C N 
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to which the phase is <J> N , and 

applying a defuzzification procedure using Y N at phases $ N and weights W N to determine the instantaneous phase 

5 [0015] Preferably, the identifiable features include zero crossings, peaks, inflection points or plateaus of the proto- 
type fiow-vs-time waveform. Furthermore, said weights can be unity, or chosen to reflect the anticipated reliability of 
deduction of the particular feature. 

[0016] The invention further discloses a method for calculating instantaneous phase in the respiratory cycle as a 
continuous variable, as described above, in which the step of calculating respiratory airflow includes a low pass filtering 
10 . step to reduce non-respiratory noise, in which the time constant of the low pass filter is an increasing function of an esti- 
mate of the length of the respiratory cycle. 

[0017] The invention further discloses a method for measuring the instantaneous phase in the respiratory cycle as 
a continuous variable as described above, in which the defuzzification step includes a correction for any phase delay 
introduced in the step of low pass filtering respiratory airflow. 
15 [0018] The invention further discloses a method for measuring the average respiratory rate, comprising the steps 
of: 

measuring leak-corrected respiratory airflow, 

from the respiratory airflow, calculating the instantaneous phase $ in the respiratory cycle as a continuous variable 
20 from 0 to 1 revolution, calculating the instantaneous rate of change of phase oWdt, and 

calculating the average respiratory rate by low pass filtering said instantaneous rate of change of phase d^dt. 

[0019] Preferably, the instantaneous phase is calculated by the methods described above. 

[0020] The invention further discloses a method for providing ventilatory assistance in a spontaneously breathing 
25 subject, comprising the steps, performed at repeated sampling intervals, of: 

ascribing a desired waveform template function n((>), with domain 0 to 1 revolution and range 0 to 1 f 
calculating the instantaneous phase $ in the respiratory cycle as a continuous variable from 0 to 1 revolution, 
selecting a desired pressure modulation amplitude A, 
30 calculating a desired instantaneous delivery pressure as an end expiratory pressure plus the desired pressure 

modulation amplitude A multiplied by the value of the waveform template function n(t» at the said calculated phase 
<t>, and 

setting delivered pressure to subject to the desired delivery pressure. 

35 [0021] The invention further discloses a method for providing ventilatory assistance in a spontaneously breathing 
subject as described above, in which the step of selecting a desired pressure modulation amplitude is a fixed amplitude. 
[0022] The invention further discloses a method for providing ventilatory assistance in a spontaneously breathing 
subject as described above, in which the step of selecting a desired pressure modulation amplitude in which said ampli- 
tude is equal to an elastance multiplied by an estimate of the subject's tidal volume. 

40 [0023] The invention further discloses a method for providing ventilatory assistance in a spontaneously breathing 
subject as described above, in which the step of selecting a desired pressure modulation amplitude comprises the sub- 
steps of: 

specifying a typical respiratory rate giving a typical cycle time, 
45 specifying a preset pressure modulation amplitude to apply at said typical respiratory rate, 

calculating the observed respiratory rate giving an observed cycle time, and 

calculating the desired amplitude of pressure modulation as said preset pressure modulation amplitude multiplied 
by said observed cycle time divided by the said specified cycle time. 

50 [0024] The invention further discloses a method for providing ventilatory assistance in a spontaneously breathing 
subject, including at least the step of determining the extent that the subject is adequately ventilated, to said extent the 
phase in the respiratory cycle is determined from the subject's respiratory airflow, but to the extent that the subject's 
ventilation is inadequate, the phase in the respiratory cycle is assumed to increase at a pre-set rate, and setting mask 
pressure as a function of said phase. 

55 [0025] The invention further discloses a method for providing ventilatory assistance in a spontaneously breathing 
subject, comprising the steps of: measuring respiratory airflow, determining the extent to which the instantaneous 
phase in the respiratory cycle can be determined from said airflow, to said extent determining said phase from said air- 
flow but to the extent that the phase in the respiratory cycle cannot be accurately determined, the phase is assumed to 
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increase at a preset rate, and delivering pressure as a function of said phase. 

[0026] The invention further discloses a method for calculating the instantaneous inspired volume of a subject, 
operable substantially without run-away under conditions of suddenly changing leak, the method comprising the steps 
of: 

5 

determining respiratory airflow approximately corrected for leak, 

calculating an index J varying from 0 to 1 equal to the fuzzy extent to which said corrected respiratory airflow is 
large positive for ledger than expected, or large negative for longer than ex 
identifying the start of inspiration, and 
10 calculating the instantaneous inspired volume as the integral of said corrected respiratory airflow multiplied by the 

fuzzy negation of said index J with respect to time, from start of inspiration. 

[0027] The invention further discloses a method "A" for providing ventilatory assistance in a spontaneously breath- 
ing subject, the method comprising the steps, performed at repeated sampling intervals, of. 

15 

determining respiratory airflow approximately corrected for leak, 

calculating an index J varying from 0 to 1 equal to the fuzzy extent to which said respiratory airflow is large positive 
for longer than expected, or large negative for longer than expected, 

calculating a modified airflow equal to said respiratory airflow multiplied by the fuzzy negation of said index J, 
20 identifying the phase in the respiratory cycle, 

calculating the instantaneous inspired volume as the integral of said modified airflow with respect to time, with the 
integral held at zero during the expiratory portion of the respiratory cycle, 

calculating a desired instantaneous delivery pressure as a function at least of the said instantaneous inspired vol- 
ume, and 

25 setting delivered pressure to subject to the desired delivery pressure. 

[0028] The invention further discloses a method H B W for providing ventilatory assistance in a spontaneously breath- 
ing subject comprising the steps of: 

30 determining respiratory airflow approximately corrected for leak, 

calculating an index J varying from 0 to 1 equal to the fuzzy extent to which the respiratory airflow is large positive 
for longer than expected, or large negative for longer than expected, 
identifying the phase in the respiratory cycle, 

calculating a modified respiratory airflow equal to the respiratory airflow multiplied by the fuzzy negation of said 

35 index J, . 

calculating the instantaneous inspired volume as the integral of the modified airflow with respect to time, with the 
integral held at zero during the expiratory portion of the respiratory cycle, 

calculating the desired instantaneous delivery pressure as an expiratory pressure plus a resistance multiplied by 
the instantaneous respiratory airflow plus a nonlinear resistance multiplied by the respiratory airflow multiplied by 
40 the absolute value of the respiratory airflow plus an elastance multiplied by the said adjusted instantaneous 
inspired volume, and 

setting delivered pressure to subject to the desired delivery pressure. 

[0029] The invention yet further discloses a method "C" for providing assisted ventilation to match the subject's 
45 need, comprising the steps of: 

describing a desired waveform template function n(<|>). with domain 0 to 1 revolution and range 0 to 1 . 
determining respiratory airflow approximately corrected for leak, 

calculating an index J varying from 0 to 1 equal to the fuzzy extent to which the respiratory airflow is targe positive 
50 for longer than expected, or large negative for longer than expected, 
calculating J peak equal to the recent peak of the index J, 
calculating the instantaneous phase in the respiratory cycle, 

calculating a desired amplitude of pressure modulation, chosen to servo-control the degree of ventilation to at least 
exceed a specified ventilation, 

55 calculating a desired delivery pressure as an end expiratory pressure plus the calculated pressure modulation 
amplitude A multiplied by the value of the waveform template function n«>) at the said calculated phase <)>. and 
setting delivered pressure to subject to said desired instantaneous delivered pressure. 
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[0030] The invention yet further discloses a method for providing assisted ventilation to match the subject's need, 
as described above, in which the step of calculating a desired amplitude of pressure modulation, chosen to servo-con- 
trol the degree of ventilation to at least exceed a specified ventilation, comprises the steps of: 

calculating a target airflow equal to twice the target ventilation divided by the target respiratory rate, 
deriving an error term equal to the absolute value of the instantaneous low pass filtered respiratory airflow minus 
the target airflow, and 

calculating the amplitude of pressure modulation as the integral of the error term multiplied by a gain, with the inte- 
gral clipped to lie between zero and a maximum. 

[0031] The invention yet further discloses a method for providing assisted ventilation to match the subject's need, 
as described above, in which the step of calculating a desired amplitude of pressure modulation, chosen to servo-con- 
trol the degree of ventilation to at least exceed a specified ventilation, comprises the following steps: 

15 calculating a target airflow equal to twice the target ventilation divided by the target respiratory rate, 

deriving an error term equal to the absolute value of the instantaneous low pass filtered respiratory airflow minus 
the target airflow, 

calculating an uncorrected amplitude of pressure modulation as the integral of the error ten multiplied by a gain, 
with the integral clipped to lie between zero and a maximum, 
20 calculating the recent average of said amplitude as the low pass filtered amplitude, with a time constant of several 

times the length of a respiratory cycle, and 

setting the actual amplitude of pressure modulation to equal the said low pass filtered amplitude multiplied by the 
recent peak jamming index J peak P |us the uncorrected amplitude multiplied by the fuzzy negation of Jpeak- 

25 [0032] The invention yet further discloses a method for providing assisted ventilation to match the subject's need, 
and with particular application to subjects with varying respiratory mechanics, insufficient respiratory drive, abnormal 
chemoreceptor reflexes, hypoventilation syndromes, or Cheyne Stokes breathing, combined with the advantages of 
proportional assist ventilation adjusted for sudden changes in leak, comprising the steps, performed at repeated sam- 
pling intervals, of: 



30 



35 



calculating the instantaneous mask pressure as described for methods "A" or "B" above, 
calculating the instantaneous mask pressure as described for method "C" above, 
calculating a weighted average of the above two pressures , and 
setting the mask pressure to the said weighted average. 



[0033] The invention yet further discloses apparatus to give effect to each one of the methods defined, including 
one or more transducers to measure flow and/or pressure, processor means to perform calculations and procedures, 
flow generators for the supply of breathable gas at a pressure above atmospheric pressure and gas delivery means to 
deliver the breathable gas to a subject's airways. 
40 [0034] The apparatus can include ventilators, ventilatory assist devices, and CPAP devices including constant level, 
bi-level or autosetting level devices. 

[0035] It is to be understood that while the algorithms embodying the invention are explained in terms of fuzzy logic, 
approximations to these algorithms can be constructed without the use of the fuzzy logic formalism. 



45 Brief Description of the Drawings 

[0036] A number of embodiments will now be described with reference to the accompanying drawings in which: 

Figs. 1a and 1b show apparatus for first and second embodiments of the invention respectively; 

50 Fig. 2 is a pressure waveform function n(<l>) used in the calculation of the desired instantaneous delivery pressure 

as a function of the instantaneous phase 0> in the respiratory cycle for a first embodiment of the invention; 
Fig 3 shows fuzzy membership functions for calculating the degree of membership in each of five magnitude fuzzy 
sets ("large negative", "small negative", "zero", "small positive", and "large positive") from the normalized respira- 
tory airflow according to the first embodiment of the invention; and 

55 Fig. 4 shows fuzzy membership functions for calculating the degree of membership in each of five rate of change 
fuzzy sets ("rising fast", "rising slowly", "steady", "falling slowly", and "falling fast") from the normalized rate of 
change of airflow according to the first embodiment of the invention; 

Fig. 5 is a pressure waveform function n(4>) used in the calculation of the desired instantaneous delivery pressure 
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as a function of the instantaneous phase <J> in the respiratory cycle for a second embodiment of the invention; 
Fig. 6 shows calculation of a quantity "lead-in" as a function of time since the most recent mask off-on transition: 
Fig. 7 shows a fuzzy membership function for fuzzy set A, as a function of time since the most recent expiratory-to- 
inspiratory (negative-to-positive) zero crossing of the respiratory airflow signal, such that the membership function 
measures the extent to which the respiratory airflow has been positive for longer than expected; 
Fig. 8 shows a membership function for fuzzy set B| as a function of respiratory airflow, such that the membership 
function measures the extent to which respiratory airflow is large positive: . 

Fig. 9 shows an electrical analog of the calculation of a recent peak jamming index J PEA K from the instantaneous 
jamming index J; 

Fig. 1 0 shows the calculation of the time constant t used in low pass filtering steps in the calculation of the conduct- 
ance of a leak, as a function of the recent peak jamming index Jpeak* 

Fig. 1 1 shows a prototypical respiratory flow-time curve, with time on the x-axis, marking nine features; 

Fig. 12 shows membership functions for fuzzy sets "large negative", "small negative", "zero", "small positive", and 

"large positive" as functions of normalized respiratory airflow according to a second embodiment of the invention; 

Fig. 13 shows membership functions for fuzzy sets 'falling", "steady", and "rising" as functions of normalized rate 

of change of respiratory airflow df/dt according to a second embodiment of the invention; 

Fig. 14 shows the membership function for fuzzy set "hypopnea"; 

Fig. 15 shows the calculation of the time constant x for calculation of normalized recent ventilation, as a function of 
"servo gain" being the gain used for servo-control of minute ventilation to at least exceed a specified target ventila- 
tion; 

Fig 1 6 shows the membership function for fuzzy set "hyperpnea" as a function of normalized recent ventilation; 
Fig 1 7 shows the membership function for fuzzy set "big leak" as a function of leak; 

Fig. 1 8 shows the membership functions for fuzzy sets "switch negative" and "switch positive" as a function of nom- 
alized respiratory airflow; 

Fig. 1 9 shows the membership functions for fuzzy sets "insp.phase" and "exp_phase" as functions of the instanta- 
neous phase in the respiratory cycle 

Fig. 20 shows schematically how function W(y), used in defuzzification, calculates the area (shaded) of an isosce- 
les triangle of unit base and height cut off below height y; 

Figs. 21-26 show actual 60 second flow and pressure tracings from the second embodiment of the invention during 
operation; the vertical scale for flow (heavy trace) is ±1 Usee, inspiration upwards and the vertical scale for the 
pressure (light trace) is 0-25 cmH 2 0; where: 

Fig. 21 shows that a short central apnea (b) is permitted when effort ceases at point (c) after a preceding deep 
breath (a); 1 

Fig. 22 shows chat a central apnea is not permitted when effort ceases at arrow (a) without a preceeding deep 
breath; 

Fig. 23 is recorded with servo gain set high, and shows that a central apnea is no longer permitted when effort 
ceases at arrow (a) despite preceding deep breathing: 

Fig. 24 shows automatically increasing end-inspiratory pressure as the subject makes voluntarily deeper inspira- 
tory efforts; 

Fig. 25 is recorded with a somewhat more square waveform selected, and shows automatically increasing pressure 
support when the subject voluntarily attempts to resist by stiffening the chest wall at point (a); 
Fig. 26 shows that with sudden onset of a sever 1 .4 Usee leak at (a), the flow signal returns to baseline (b) within 
the span of a single breath, and pressure continues to cycle correctly throughout; and 

Fig. 27 shows an actual 60 second tracing showing respiratory airflow (heavy trace, ± 1 L/sec full scale) and instan- 
taneous phase (light trace, 0-1 revolution full scale). 

Description of Preferred Embodiments 

[0037] The two embodiments to be described are ventilators that operate in a manner that seeks to simultaneously 
achieve the three goals stated above. 



First Embodiment 

[0038] Apparatus to give effect to a first embodiment of the apparatus is shown in Fig. 1a. A blower 10 supplies a 
breathable gas to mask 11 in communication with the subject's airway via a delivery tube 12 and exhausted via a 
exhaust diffuser 1 3. Airflow to the mask 1 1 is measured using a pneumotachograph 1 4 and a differential pressure trans- 
ducer 15. The mask flow signal from the transducer 15 is then sampled by a microprocessor 16. Mask pressure is 
measured at the port 17 using a pressure transducer 18. The pressure signal from the transducer 18 is then sampled 
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by the microprocessor 16. The microprocessor 16 sends an instantaneous mask pressure request signal to the servo 
19, which compares said pressure request signal with actual pressure signal from the transducer 18 to the control fan 
motor 20. The microprocessor settings can be adjusted via a serial port 21 . 

[0039] It is to be understood tat the mask could equally be replaced with a tracheotomy tube, endotracheal tube, 
nasal pillows, or other means of making a sealed connection between the air delivery means and the subject's airway. 
[0040] The microprocessor 1 6 Is programmed to perform the following steps, to be considered in conjunction with 
Tables 1 and 2. 



Table 1 





Fuzzy Inference Rules for a first embodiment 




N 


Fuzzy Interference Rule 


Fuzzy Phase 


15 


1 
2 


if size is 
if size is 


Zero 

Small Positive 


and rate of 

and rate of change 
is 


Increasing 
Increasing Slowly 


then phase is 
then phase is 


Start Inspiration 
Early Inspiration 




\j 


if size is 


Large Positive 


and rate of change 
is 


Steady 


then phase is 


Peak Inspiration 


20 


4 


if size is 


Small Positive 


and rate of change 
is 


Decreasing Slowly 


then phase is 


Late Inspiration 




5 


if size is 


Zero 


and rate of change 
is 


Decreasing Fast 


then phase is 


Start Expiration 


25 


6 


if size is 


Small Negative 


and rate of change 
is 


Decreasing Slowly 


then phase is 


Early Expiration 




7 


if size is 


Large Negative 


and rate of change 
is 


Steady 


then phase is 


Peak Expiration 


30 


8 


if size is 


Small Negative 


and rate of change 
is 


Increasing Slowly 


then phase is 


Late Expiration 




9 


if size is 


Zero 


and rate of change 
is 


Steady 


then phase is 


Expiratory Pause 


35 


10 


always 






phase is 


Unchanged 



Table 2 



Association of phases with fuzzy rules for a first embodiment 


N 


Phase 




1 


Start Inspiration 


0.0 


2 


Early Inspiration 


values intermediate between 0.0 and 0.5 


3 


Peak Inspiration 


4 


Late Inspiration 


5 


Start Expiration 


0.50 


6 


Early Expiration 


values intermediate between 0.5 and 1.0 


7 


Peak Expiration 


8 


Late Expiration 


9 


Expiratory Pause 
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Table 2 (continued) 



Association of phases with fuzzy rules for a first embodiment. 


N 


Phase 




10 


Unchanged 


0 



10 1 . Set desired target values for the duration of inspiration TI TGT , duration of expiration TEtqt, and minute ventilation 

V TGT . Choose suitable constants P 0 and A STD where P 0 is the desired end expiratory pressure, and A STD is the 
desired increase in pressure above P 0 at end inspiration for a breath of duration TT TGT =TI TGT +TE TGT . 

2. Choose a suitable pressure waveform function n(4>), such as that shown in Fig. 2, such that the desired delivery 
15 pressure at phase 4> will be given by: 

p = p 0 + An(<D) 

where the amplitude A equals the difference between the end inspiratory pressure and end expiratory pressure. 
20 However, other waveforms may be suitable for subjects with particular needs. 

3. Initialize the phase O in the respiratory cycle to zero, and initialize the current estimates of actual inspiratory and 
expiratory duration Tl and TE to TI TGT and TE TGT respectively. 

25 4. Initialize the rate of change of phase during inspiration A4>| between sampling intervals of length T to: 

A4>+ = 0.5 T / Tl TGT 

5. Initialize the rate of change of phase during expiration A<J>£ to: 



30 



AO E = 0.5T/TE TGT 



6. Measure the instantaneous respiratory airflow f RESP 
35 7. Calculate the average total breath duration TT = Tl + TE 

8. Low pass filter the respiratory airflow with an adjustable time constant xf, where xf is a fixed small fraction of TT. 

9. Calculate the instantaneous ventilation V, as half the absolute value of the respiratory airflow: 

40 

V = 0.5 |f RESP | 

10. From the target ventilation V TGT and the measured minute ventilation V, derive an error term V ERRt such that 
large values of V ERR indicate inadequate ventilation: 

V ERR = /(V TGr V)dt 

1 1 . Take V BAR as the result of low pass filtering V with a time constant iVbar which is long compared with TT. 
so 12. Calculate a normalized airflow f|sjoRM» where 

^ NORM = ^RESP^BAR* 

13. From f^oRM. calculate the degree of membership in each of the fuzzy sets whose membership functions are 
55 shown in Fig. 3. 

14. Calculate a normalized rate of change df NO RM/ d< & equal to df NO RM/dt divided by the current estimate of the 
average respiratory cycle time TT. 
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1 5. From the normalized rate of change, calculate the degree of membership In each of the fuzzy sets shown In Fig. 



1 6. For each row N in Table 1 f calculate the degree of membership g N in the fuzzy set shown in the column labelled 
5 Fuzzy Phase, by applying the fuzzy inference rules shown. 

1 7. Associate with the result of each of the N rules a phase On as shown in Table 2, noting that O 10 is the current 
phase O 

10 18. Increase each of the On excepting O 10 by 0.89 x/TT, to compensate for the previous low pass filtering step. 

1 9. Calculate a new instantaneous phase O IN st as the angle to the center of gravity of N unit masses at polar coor- 
dinates of radius g^ and angle On revolutions. 

15 20. Calculate the smallest signed difference AO| NST bewteen the phase estimated in the previous step and the cur- 

rent phase. 

A4> inst = 1 - (AO lNST -O) (O )NST - 0 > 0.5) 

20 ^ INST = <*> INST - * + 1 INST " * < ' 0 5 ) 

AOINST = O INST - <D (otherwise) 

21 . Derive a revised estimate AO REV equal to a weighted mean of the value calculated in the previous step and the 
25 average value (AO ( or AO E as appropriate). 

A4> = (1-W) AO , + WAO |NST (0 < O < 0.5) 

A<J> = (1-W) A* | + WAO inst (otherwise) 

30 

Smaller values of W will cause better tracking of phase if the subject is breathing regularly, and larger values will 
cause better tracking of phase if the subject is breathing irregularly. 

22. Derive a blending fraction B, such that the blending fraction is unity if the subject's ventilation is well above V TGTt 
35 zero if the subject is breathing near or below V TGT , and increasing proportionally from zero to unity as the subject's 

ventilation increases through an intermediate range. 

23. Calculate AO BLEN d influenced chiefly by AO calculated in step 21 from the subject's respiratory activity rf the 
subject's ventilation is well above V TGT ; influenced chiefly by the target respiratory duration if the subject is breath- 

40 ing near or below V TGT ; and proportionally between these two amounts if ventilation is in an intermediate range: 

AO BLEND = B AO+ °- 5 ( 1 B ) T ' TI TGT (0 < O < 0.5) 

A4> BLEND = B AO+ 0.5 (1 -B) T / TE TGT (otherwise) 

45 

24. Increment 0 by A<|>bl E nd 

25. Update the avenge rate of change of phase (AO) or AO E as appropriate). 

50 A^I=T/t vbar (AO blend -AO,) (0 < O < 0.5) 

AO E = T/x vbar ( A<p blend " A<t> e> (otherwise) 

26. Recalculate the approximate duration of inspiration TI and expiration TE: 

55 

TI = 0.5T/ AO, 
TE = 0.5T/AO E 



10 



10 



EP 1 005 829 A1 

27. Calculate the desired mask pressure modulation amplitude Aq: 

A D = A STD /2 (TT<TT STD /2) 
a d = 2 * a std (TT>2 • TT STD ) 

a d = a std • tt/tt std (otherwise). 

28. From the error term V ERR , calculate an additional mask pressure modulation amplitude Ag: 

A E = K*V ERR (forV ERR >0) 
A E = 0 (otherwise) 

is where larger values of K will produce a faster but less stable control of the degree of assistance, and smaller values 

of K will produce slower but more stable control of the degree of assistance. 

29. Set the mask pressure Pmask to: 

20 Pmask = Po + ( A D +A E )n(<l>) 

30. Wait for a sampling interval T t short compared with the duration of a respiratory cycle, and then continue at the 
step of measuring respiratory air flow. 

25 Measurement of respiratory airflow 



[0041] As follows from above, it is necessary to respiratory airflow, which is a standard procedure to one skilled in 
the art. In the absence of leak, respiratory airflow can be measured directly with a pneumotachograph placed between 
the mask and the exhaust. In the presence of a possible leak, one method disclosed in European Publication No 0 651 
30 971 incorporated herein by cross-reference is to calculate the mean flow through the leak, and thence calculate the 
amount of modulation of the pneumotachograph flow signal due to modulation of the flow through the leak induced by 
changing mask pressure, using the following steps: 

1 . Measure the airflow at the mask f^ASK using a pneumotachograph 
35 2. Measure the pressure at the mask Pmask 

3. Calculate the mean leak as the low pass filtered airflow, with a time constant long compared with a breath. 

4. Calculate the mean mask pressure as the low pass filtered mask pressure, with a time constant long compared 
with a breath. 

5. Calculate the modulation of the flow through the leak as: 

40 

- S(leak) = 0.5 times the mean leak times the inducing pressure, 
where the inducing pressure is Pmask ~ mean mask pressure. 



45 Thence the instantaneous respiratory airflow can be calculated as: 



f resp = f mask * mean ,eak - 5(leak) 

A convenient extension as further disclosed in EP 0 651 971 (incorporated herein by cross-reference) is to measure air- 
so flow f turbine ana " pressure Pturbine at ^ e outlet of the turbine, and thence calculate Pmask and Wsk bv allowing for 
the pressure drop down the air delivery hose, and the airflow lost via the exhaust: 

1 - AP HOS E = K l( F TURBtNE)" K 2( F TURBINE) 

2 - PMASK = P TURBINE" AP HOSE 
55 3 * F EXHAUST = K3 ^ P MASK 

4* F MASK = ^ TURBINE " ^ EXHAUST 
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Alternative embodiment 

[0042] The following embodiment is particularly applicable to subjects with varying respiratory mechanics, insuffi- 
cient respiratory drive, abnormal chemoreceptor reflexes, hypoventilation syndromes, or Cheyne Stokes breathing, or 

5 to subjects with abnormalities of the upper or lower airways, lungs, chest wall, or neuromuscular system. 

[0043] Many patients with severe lung disease cannot easily be treated using a smooth physiological pressure 
waveform, because the peak pressure required is unacceptably high, or unachievable with for example a nose-mask. 
Such patients may prefer a square pressure waveform, in which pressure rises explosively fast at the moment of com- 
mencement of inspiratory effort This may be particularly important in patients with high intrinsic PEER in which it is not 

10 practicable to overcome the intrinsic PEEP by the use of high levels of extrinsic PEEP or CPAP, due to the risk of hyper- 
inflation. In such subjects, any delay in triggering is perceived as very distressing, because of the enormous mis-match 
between expected and observed support. Smooth waveforms exaggerate the perceived delay, because of the time 
taken for the administered pressure to exceed the intrinsic PEEP. This embodiment permits the use of waveforms var- 
ying continuously from square (suitable for patients with for example severe lung or chest wall disease or high intrinsic 

15 PEEP) to very smooth, suitable for patients with normal lungs and chest wall, but abnormal respiratory control, or neu- 
romuscular abnormalities. This waveform is combined either with or without elements of proportional assist ventilation 
(corrected for sudden changes in leak), with servo-control of the minute ventilation to equal or exceed a target ventila- 
tion. The latter servo-control has an adjustable gain, so that subjects with for example Cheyne Stokes breathing can be 
treated using a very high servo gain to over-ride their own waxing and waning patterns; subjects with various central 

20 hypoventilation syndromes can be treated wit a low servo gain, so that short central apneas are permitted, for example 
to cough, clear the throat, talk, or roll over in bed, but only if they follow a previous period of high ventilation; and normal 
subjects are treated with an intermediate gain. 
[0044] Restating the above in other words: 

25 • The integral gain of the servo-control of the degree of assistance is adjustable from very fast (0.3 
cmH 2 0/L/sec/sec) to very slow. Patients with Cheyne-Stokes breathing have a very high ventilatory control loop 
gain, but a long control loop delay, leading to hunting. By setting the loop gain even higher, the patient's controller 
is stabilized. This prevents the extreme breathlessness that normally occurs during each cycle of Cheyne-Stokes 
breathing, and this is very reassuring to the patient. It is impossible for them to have a central apnea. Conversely, 

30 subjects with obesity-hypoventilation syndrome have low or zero loop gain. They will not feel breathless during a 

central apnea. However, they have much mucus and need to cough, and are also often very fidgety, needing to roll 
about in bed. This requires that they have central apneas which the machine does not attempt to treat. By setting 
the loop gain very low, the patient is permitted to take a couple of deep breaths and then have a moderate-length 
central apnea while coughing, rolling over, etc, but prolonged sustained apneas or hypopneas are prevented. 

35 

Sudden changes in leakage flow are detected and handled using a fuzzy logic algorithm. The principle of the algo- 
rithm is that the leak filter time constant is reduced dynamically to the fuzzy extent that the apparent respiratory air- 
flow is a long way from zero for a long time compared with the patient's expected respiratory cycle length. 

40 • Rather than simply triggering between two states (IPAP, EPAP), the device uses a fuzzy logic algorithm to estimate 
the position in the respiratory cycle as a continuous variable. The algorithm permits the smooth pressure waveform 
to adjust it's rise time automatically to the patient's instantaneous respiratory pattern. 

The fuzzy phase detection algorithm under normal conditions closely tracks the patient's breathing. To the extent 
45 that there is a high or suddenly changing leak, or the patient's ventilation is low, the rate of change of phase (res- 

piratory rate) smoothly reverts to the specified target respiratory rate. Longer or deeper hypopneas are permitted 
to the extent that ventilation is on average adequate. To the extent that the servo gain is set high to prevent Cheyne 
Stokes breathing, shorter and shallower pauses are permitted. 

so • Airflow filtering uses an adaptive filter, which shortens it's time constant if the subject is breathing rapidly, to give 
very fast response times, and lenthens if the subject is breathing slowly, to help eliminate cardiogenic artifact. 

The fuzzy changing leak detection algorithm, the fuzzy phase detection algorithm with its differential handling of 
brief expiratory pauses, and handling of changing leak, together with the smooth waveform severally and coopera- 
55 tively make the system relatively immune to the effects of sudden leaks. 

• By suitably setting various parameters, the system can operate in CPAP, bilevel spontaneous, bilevel timed, propor- 
tional assist ventilation, volume cycled ventilation, and volume cycled servo-ventilation, and therefore all these 
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modes are subsets of the present embodiment However, the present embodiment permits states of operation that 
can not be achieved by any of the above states, and is therefore distinct from them. 

Notes 

5 

[0045] 

Note 1: in this second embodiment, the names and symbols used for various quantities may be different to those 
used in the first embodiment 

w Note 2: The term "swing" is used to refer to the difference between desired instantaneous pressure at end inspira- 

tion and the desired instantaneous pressure at end expiration. 

Note 3: A fuzzy membership function is taken as returning a value between zero for complete nonmembership and 
unity for complete membership. Fuzzy intersection A AND B is the lesser of A and B, fuzzy union A OR B is the 
larger of A and B, and fuzzy negation NOT A is 1 -A. 
15 Note 4: root(x) is the square root of x, abs(x) is the absolute value ofx t sign(x) is - 1 ifx is negative, and + 1 other- 

wise. An asterisk (*) is used to explicitly indicate multiplication where this might not be obvious from context. 

Apparatus 

20 [0046] The apparatus for the second embodiment is shown in Fig. 1b. The blower 110 delivers air under pressure 
to the mask 11 1 via the air delivery hose 1 12. Exhaled air is exhausted via the exhaust 1 1 3 in the mask 1 1 1 . The pneu- 
motachograph 114 and a differential pressure transducer 115 measure the airflow in the nose 112. The flow signal is 
delivered to the microprocessor 116. Pressure at any convenient point 117 along the nose 112 is measured using a 
pressure transducer 118. The output from the pressure transducer 118 is delivered to the microcontroller 1 16 and also 

25 to a motor servo 119. The microprocessor 116 supplies the motor servo 119 with a pressure request signal, which is 
then compared with the signal from the pressure transducer 118 to control the blower motor 120. User configurable 
parameters are loaded into the microprocessor 1 16 via a communications port 121 , and the computed mask pressure 
and flow can if desired be output via the communications port 121. 

30 Initialization 

[0047] The following user adjustable parameters are specified and stored: 



max permissible pressure 


maximum permissible mask pressure 


max swing 


maximum permissible difference between end inspiratory pressure and end 
expiratory pressure. 


min swing 


minimum permissible difference between end inspiratory pressure and end 
expiratory pressure. 


epap 


end expiratory pressure 


min permissible pressure 


minimum permissible mask pressure 


target ventilation 


minute ventilation is sevo-controlled to equal or exceed this quantity 


target frequency 


Expected respiratory rate. If the patient is achieving no respiratory airflow, 
the pressure will cycle at this frequency. 


target duty cycle 


Expected ratio of inspiratory time to cycle time. If the patient is achieving no 
respiratory airflow, the pressure will foPow this duty cycle. 


linear resistance andquad resistance 


resistive unloading ^linear resistance * f + quad.resistance * f 2 
sign(f),where f is the respiratory airflow.where sign(x) = -1 for x < 0, +1 oth- 
erwise 


elastance 


Unload at least this much elastance 


servo gain 


gain for servo-control of minute ventilation to at least exceed target ventila- 
tion. 
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(continued) 



max permissible pressure 


maximum permissible mask pressure 


waveform time constant 


Elastic unloading waveform time constant as a fraction of inspiratory dura- 
tion. (0.0 = square wave) 


hose resistance 


AP from pressure sensing port to inside mask = hose resistance times the 
square of the flow in the Intervening tubing. 


diffuser conductance 


Flow through the mask exhaust port = diffuser conductance * root mask 
pressure 



10 
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[0048] 
[0049] 



[0050] 

tively to: 



At initialization, the following are calculated from the above user-specified settings: 

The expected duration of a respiratory cycle, of an inspiration, and of an expiration are set respectively to: 
STD T TOT - 60 I target respiratory rate 
STD T, = STD T 

jot * ^'P©* o*ury cycle 
STD T E = STD T TOT - STD 7", 
The standard rates of change of phase (revolutions per sec) during inspiration and expiration are set respec- 



STD cty , = 0.5 1 STD T t 

STD o*4> 5 = 0.51 STD T E 
[0051] The instantaneous elastic support at any phase <|> in the respiratory cycle is given by: 

PEL(<(>) = swing * n(<|>) 
where swing is the pressure at end inspiration minus the pressure at end expiration. 



n(<|») = e t<(> 
e^-O.S) 



during inspiration, 
during expiration 



and t is the user-selectable waveform time constant. 

If x - 0, then n(<|>) is a square wave. The maximum implemented value for t = 0.3, producing a waveform approximately 
40 as shown in Fig. 5. 

[0052] The mean value of n($) is calculated as follows: 



45 



Jn(4>)cty 



50 

Operations Performed every 20 Milliseconds 

[0053] The following is an overview of routine processing done at 50 Hz: 

55 measure flow at flow sensor and pressure at pressure sensing port 

calculate mask pressure and flow from sensor pressure and flow 
calculate conductance of mask leak 
calculate instantaneous airflow through leak 
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calculate respiratory airflow and low pass filtered respiratory airflow 
calculate mask on-off status and lead-in 
calculate instantaneous and recent peak jamming 
calculate time constant for leak conductance calculations 
5 calculate phase in respiratory cycle 

update mean rates of change of phase for inspiration and expiration, lengths of inspiratory and expiratory times, 
and respiratory rate 

add hose pressure loss to EPAP pressure 
add resistive unloading 
10 calculate instantaneous elastic assistance required to servo-control ventilation 

estimate instantaneous elastic recoil pressure using various assumptions 
weight and combine estimates 
add servo pressure to yield desired sensor pressure 
servo-control motor speed to achieve desired sensor pressure 

15 

[0054] The details of each step will now be explained. 

Measurement of Flow and Pressure 

20 [0055] Flow is measured at the outlet of the blower using a pneumotachograph and differential pressure transducer. 
Pressure is measured at any convenient point between the blower outlet and the mask. A humidifier and/or anti-bacte- 
rial filter may be inserted between the pressure sensing port and the blower. Flow and pressure are digitized at 50 Hz 
using an A/D converter. 

25 Calculation of mask flow and pressure 

[0056] The pressure loss from pressure measuring point to mask is calculated from the flow at the blower and the 
(quadratic) resistance from measuring point to mask. 

30 Hose pressure loss = sign(flow) * hose resistance * flow 

where sign(x) = -1 for x < 0, +1 otherwise. The mask pressure is then calculated by subtracting the hose pressure loss 
from the measured sensor pressure: 

35 Mask pressure - sensor pressure - hose pressure loss 

[0057] The flow through the mask exhaust diffuser is calculated from the known parabolic resistance of the diffuser 
holes, and the square root of the mask pressure: 

40 diffuser flow - exhaust resistance * signfmask pressure) * root(abs(mask pressure)) 

[0058] Finally, the mask flow is calculated: 

mask flow = sensor flow - diffuser flow 

45 

[0059] The foregoing describes calculation of mask pressure and flow in the various treatment modes. In diagnostic 
mode, the patient is wearing only nasal cannulae, not a mask. The cannula is plugged into the pressure sensing port 
The nasal airflow is calculated from the pressure, after a linearization step, and the mask pressure is set to zero by def- 
inition. 

so 

Conductance of leak 

[0060] The conductance of the leak is calculated as follows: 
55 root mask pressure = sign (P ^sk) V abs ( p mask) 

LP mask airflow = low pass filtered mask airflow 
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LP root mask pressure = low pass filtered root mask pressure 
conductance of leak = LP mask airflow / LP root mask pressure 

5 [0061] The time constant for the two low pass filtering steps is initialized to 10 seconds and adjusted dynamically 
thereafter (see below). 

Instantaneous flow through leak 

10 [0062] The instantaneous flow through the leak is calculated from the instantaneous mask pressure and the con- 
ductance of the leak: 

instantaneous leak = conductance of leak * root mask pressure 

15 Respiratory Airflow 

[0063] The respiratory airflow is the difference between the flow at the mask and the instantaneous leak: 

respiratory airflow = mask flow - instantaneous leak 

20 

Low pass Altered respiratory airflow 

[0064] Low pass filter the respiratory airflow to remove cardiogenic airflow and other noise. The time constant is 
dynamically adjusted to be 1/40 of the current estimated length of the respiratory cycle T TOT (initialized to STD_T TO t 
25 and updated below). This means that at high respiratory rates, there is only a short phase delay introduced by the filter, 
but at low respiratory rates, there is good rejection of cardiogenic airflow. 

Mask on/off status 

30 [0065] The mask is assumed to initially be off. An off-on transition is taken as occurring when the respiratory airflow 
first goes above 0.2 Usee, and an on-off transition is taken as occurring if the mask pressure is less than 2 cmH 2 0 for 
more than 1 .5 seconds. 

Lead-in 

35 ' 

[0066] Lead-in is a quantity that runs from zero if the mask is off, or has just been donned, to 1 .0 if the mask has 
been on for 20 seconds or more, as shown in Figure 6. 

Calculation of instantaneous jamming index, J 

40 

[0067] J is the fuzzy extent to which the impedance of the leak has suddenly changed. It is calculated as the fuzzy 
extent to which the absolute magnitude of the respiratory airflow is large for longer than expected. 
[0068] The fuzzy extent A| to which the airflow has been positive for longer than expected is calculated from the 
time t a since the last positive-going zero crossing of the calculated respiratory airflow signal, and the expected duration 
45 STD T| of a normal inspiration for the particular subject, using the fuzzy membership function shown in Figure 7. 

[0069] The fuzzy extent B, to which the airflow is large and positive is calculated from the instantaneous respiratory 
airflow using the fuzzy membership function shown in Figure 8. 

[0070] The fuzzy extent l ( to which the leak has suddenly increased is calculated by calculating the fuzzy intersec- 
tion (lesser) of Aj and B ( . 

so [0071] Precisely symmetrical calculations are performed for expiration, deriving l E as the fuzzy extent to which the 
leak has suddenly decreased. Ag is calculated from Tze and T E , B E is calculated from minus fpesp a nd 'e is the fuzzy 
intersection of A E and B E . The instantaneous jamming index J is calculated as the fuzzy union (larger) of indices l ( and 

55 Recent peak Jamming 

[0072] If the instantaneous jamming index is larger than the current value of the recent peak jamming index, then 
the recent peak jamming index is set to equal the instantaneous jamming index. Otherwise, the recent peak jamming 
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index is set to equal the instantaneous jamming index low pass filtered with a time constant of 1 0 seconds. An electrical 
analogy of the calculation is shown in Figure 9. 

Time constant for leak conductance calculations 

5 

[0073] If the conductance of the leak suddenly changes, then the calculated conductance will initially be Incorrect, 
and will gradually approach the correct value at a rate which will be slow if the time constant of the low pass filters is 
long, and fast if the time constant is short. Conversely, if the impedance of the leak is steady, the longer the time con- 
stant the more accurate the calculation of the instantaneous leak. Therefore, it is desirable to lengthen the time constant 

10 to the extent that the leak is steady, reduce the time constant to the extent that the leak has suddenly changed, and to 
use intermediately longer or shorter time constants if it is intermediately the case that the leak is steady. 
[0074] If there is a large and sudden increase in the conductance of the leak, then the calculated respiratory airflow 
will be incorrect. In particular, during apparent inspiration, the calculated respiratory airflow will be large positive for a 
time that is large compared with the expected duration of a normal inspiration. Conversely, if there is a sudden decrease 

15 in conductance of the leak, then during apparent expiration the calculated respiratory airflow will be large negative for 
a time that is large compared with the duration of normal expiration. 

[0075] Therefore, the time constant for the calculation of the conductance of the leak is adjusted depending on 
j peak. which is a measure of the fuzzy extent that the leak has recently suddenly changed, as shown in Figure 1 0. 
[0076] In operation, to the extent that there has recently been a sudden and large change in the leak, J peak will ^ 

20 large, and the time constant for the calculation of the conductance of the leak will be small, allowing rapid convergence 
on the new value of the leakage conductance. Conversely, if the leak is steady for a long time, J peak will be small, and 
the time constant for calculation of the leakage conductance will be large, enabling accurate calculation of the instan- 
taneous respiratory airflow. In the spectrum of intermediate situations, where the calculated instantaneous respiratory 
airflow is larger and for longer periods, J peak will be progressively larger, and the time constant for the calculation of 

25 the leak will progressively reduce. For example, at a moment in time where it is uncertain whether the leak is in fact con- 
stant, and the subject has merely commenced a large sigh, or whether in fact there has been a sudden increase in the 
leak, the index will be of an intermediate value, and the time constant for calculation of the impedance of the leak will 
also be of an intermediate value. The advantage is that some corrective action will occur very early, but without momen- 
tary total loss of knowledge of the impedance of the leak. 
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Instantaneous phase in respiratory cycle 



[0077] The current phase $ runs from 0 for start of inspiration to 0.5 for start of expiration to 1 .0 for end expiration 
= start of next inspiration. Nine separate features (peaks, zero crossings, plateaux, and some intermediate points) are 
35 identified on the waveform, as shown in Figure 1 1 . 

Calculation of normalized respiratory airflow 

[0078] The filtered respiratory airflow is normalized with respect to the user specified target ventilation as follows: 

40 

standard airflow = target ventilation / 7,5 L/min 

f = filtered respiratory airflow / standard airflow 

45 [0079] Next, the fuzzy membership in fuzzy sets large negative, small negative, zero, small positive, and large 
positive, describing the instantaneous airflow is calculated using the membership functions shown in Figure 12. For 
example, if the normalized airflow is 0.25, then the airflow is large negative to extent 0.0, small negative to extent 0.0, 
zero to extent 0.5, small positive to extent 0.5, large positive to extent 0.00. 

so Calculation of normalized rate of change of airflow 

[0080] The rate of change of filtered respiratory airflow is calculated and normalized to a target ventilation of 7.5 
L/min at 15 breaths/min as follows: 

55 standard df/dt = standard airflow * target frequency/ 15 

calculate d( filtered airfiow)/dt 
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low pass filter with a time constant of 8/50 seconds 



normalize by dividing by standard df/dt 

5 [0081] Now evaluate the membership of normalized df/dt in the fuzzy sets falling, steady, and rising, whose mem- 
bership functions are shown in Figure 13. 

Calculation of ventilation, normalized ventilation, and hypopnea 

w [0082] 

ventilation = absfrespiratory airflow) , 

low pass filtered with a time constant of STD T jqj . 

normalized ventilation = ventilation I standard airflow 

[0083] Hypopnea is the fuzzy extent to which the normalized ventilation is zero. The membership function for 
hypopnea is shown in Fig. 14. 
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Calculation of recent ventilation, normalized recent ventilation, and hyperpnea 



[0084] Recent ventilation is also a low pass filtered abs(respiratory airflow), but filtered with an adjustable time 
constant, calculated from servo gain (specified by the user) as shown in Figure 15. For example, if the servo gain is 
25 set to the maximum value of 0.3, the time constant is zero, and recent ventilation equals instantaneous abs(respira- 
tory airflow). Conversely, if servo gain is zero, the time constant is twice STD T TOT , the expected length of a typical 
breath. 



Ta rget absolute airflow - 2 * target ventilation 
normalized recent ventilation = recent ventilation / target absolute airflow 



[0085] Hyperpnea is the fuzzy extent to which the recent ventilation is large. The membership function for hyper- 
pnea is shown in Fig. 16. 

35 

Big Leak 

[0086] The fuzzy extent to which there is a big leak is calculated from the membership function shown in Figure 1 7. 

40 Additional fuzzy sets concerned with fuzzy "triggering" 

[0087] Membership in fuzzy sets switch negative and switch positive are calculated from the normalized respi- 
ratory airflow using the membership functions shown in Figure 18, and membership in fuzzy sets insp_phase and 
exp_phase are calculated from the current phase fusing the membership functions shown in Fig. 19. 

45 

Fuzzy Inference Rules for Phase 

[0088] Procedure W(y) calculates the area of an isosceles triangle of unit height and unit base, truncated at height 
y as shown in Figure 20. In the calculations that follow, recall that fuzzy intersection a AND b is the smaller of a and b, 
so fuzzy union a OR b is the larger of a and b, and fuzzy negation NOT a is 1-a. 

[0089] The first fuzzy rule indicates that lacking any other information the phase is to increase at a standard rate. 
This rule is unconditionally true, and has a very heavy weighting, especially if there is a large leak, or there has recently 
been a sudden change in the leak, or there is a hypopnea. 



55 ^STANDARD = 8+16* 

J peak + 16 * byopopnea + 16 * big leak 

[0090] The next batch of fuzzy rules correspond to the detection of various features of a typical flow-vs-time curve. 
These rules all have unit weighting, and are conditional upon the fuzzy membership in the indicated sets: 
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W early insp = W(rise and small positive) 
w peak insp - W(\arge positive AND steady AND NOT recent peak jamming) 
5 w late insp = WffctlAND small positive) 

^ early exp ~ Wffoll AND small negative) 
w peak exp = W(large negative AND steady) 

10 

w late exp = W(rise AND small negative) 

[0091] The next rule indicates that there is a legitimate expiratory pause {as opposed to an apnea) if there has been 
a recent hyperpnea and the leak has not recently changed: 

15 

w pa use - (hyperpnea AND NOT J PEAK ) * Wfsteady AND zero) 

[0092] Recalling that the time constant for hyperpnea gets shorter as servo gain increases, the permitted length 
of expiratory pause gets shorter and shorter as the servo gain increases, and becomes zero at maximum servo gain. 
The rationale for this is that (i) high servo gain plus long pauses in breathing will result in "hunting 11 of the servo-control- 
ler, and (ii) in general high servo gain is used if the subject's chemoreceptor responses are very brisk, and suppression 
of long apneas or hypopneas will help prevent the subject's own internal servo-control from hunting, thereby helping 
prevent Cheyne-Stokes breathing. 

[0093] Finally, there are two phase-switching rules. During regular quiet breathing at roughly the expected rate, 
these rules should not strongly activate, but they are there to handle irregular breathing or breathing at unusual rates. 
They have very heavy weightings. 

^ trig insp " 32 W(expiratory phase AND switch positive) 
30 w trig exp ~ 32 vVfinspiratory phase AND switch negative) 

Defuzzification 

[0094] For each of the ten fuzzy rules above, we attach phase angles fN, as shown in Table ZZZ. Note that <|>N are 
35 in revolutions, not radians. We now place the ten masses W(N) calculated above at the appropriate phase angles <|> N 
around the unit circle, and take the centroid. 



40 



45 



50 



55 



Rule 


N 




STANDARD 


1 


current § 


TRIG INSP 


2 


0.00 


EARLY INSP 


3 


0.10 


PEAK INSP 


4 


0.30 


LATE INSP 


5 


0.50 


TRIG EXP 


6 


0.5 + 0.05 k 


EARLY EXP 


7 


0.5 + 0.10 k 


PEAK EXP 


8 


0.5 + 0.20 k 


LATE EXP 


9 


0.5 + 0.4 k 


EXP PAUSE 


10 


0.5 + 0.5 k 



where k = STD T, / STD T E . 

[0095] Note that if the user has entered very short duty cycle, k will be small. For example a normal duty cycle is 
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40%, giving k = 40/60 = 0.67. Thus the expiratory peak wDI be associated with a phase angle of 0.5+0.2*0.67=0.63, cor- 
responding 26% of the way into expiratory time, and the expiratory pause would start at 0.5+0.5*0.67=0.83, corre- 
sponding to 67% of the way into expiratory time. Conversely, if the duty cycle is set to 20% in a patient with severe 
obstructive lung disease, features 6 through 10 will be skewed or compressed into earty expiration, generating an 
appropriately longer expiratory pause. 

[0096] The new estimate of the phase is the centroid, in polar coordinates, of the above ten rules: 

ZW N sin6 N 
centroid = arctan . . -^-) 

^ZWnCOS^n' 

[0097] The change in phase d<> from the current phase <j> to the centroid is calculated in polar coordinates. Thus if 
the centroid is 0.01 and the current phase is 0.99, the change in phase is d$ - 0.02. Conversely, if the centroid is 0.99 
and the current phase is 0.01 , then d<|t = -0.02. The new phase is then set to the centroid: 



4> = centroid 

[0098] This concludes the calculation of the instantaneous phase in the respiratory cycle <J>. 
20 Estimated mean duration of inspiration, expiration, cycle time, and respiratory rate 

[0099] If the current phase is inspiratory < 0.5) the estimated duration of inspiration T ( is updated: 
LP(d§ j) = low pass filtered cty with a time constant of4*STD T TOT 

25 

Clip LP(d$ t ) to the range (O.S/STD T,)/2 to 4(0.5/STD T t ) 
T,= 0.5 /clipped LP(d$ I) 

30 [0100] Conversely, if the current phase is expiratory, (<► > = 0.5) the estimated duration of expiration T E is updated: 

LP(d$ £ ) - low pass filtered d§ with a time constant of4*STD T TOT 
Clip LP(d$E) to the range (O.S/STD T E y2io 4(0.5/STDT E ) 

35 

T E = 0.5 / clipped LP(d$ E ) 

[0101] The purpose of the clipping is firstly to prevent division by zero, and also so that the calculated T ( and T E are 
never more than a factor of 4 shorter or a factor of 2 longer than expected. 
40 [0102] Finally, the observed mean duration of a breath T TOT and respiratory rate RR are: 

t tot = T, + T E 



RR = 60IT TOT 



45 

Resistive unloading 



[0103] The resistive unloading is the pressure drop across the patient's upper and lower airways, calculated from 
the respiratory airflow and resistance values stored in SRAM 

so 

f = respiratory airflow truncated to +/- 2 Usee 

resistive unloading = airway resistance * f + upper airway resistance * f 2 * sign(f) 

55 Instantaneous Elastic Assistance 

[0104] The purpose of the instantaneous elastic assistance is to provide a pressure which balances some or all of 
the elastic deflating pressure supplied by the springiness of the lungs and chest wall (instantaneous elastic pressure), 
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plus an additional component required to servo-control the minute ventilation to at least exceed on average a pre-set 
target ventilation. In addition, a minimum swing, always present, is added to the total. The user-specified parameter 
elastance is preset to say 50-75% of the known or estimated elastance of the patient's lung and chest wall. The various 
components are calculated as follows: 

5 

Instantaneous assistance based on minimum pressure swing set by physician: 
[0105] 

w instantaneous minimum assistance = minimum swing * Yl($) 

Elastic assistance required to servo-control ventilation to equal or exceed target 

[0106] The quantity servo swing is the additional pressure modulation amplitude required to servo-control the 
15 minute ventilation to at least equal on average a pre-set target ventilation. 

[0107] Minute ventilation is defined as the total number of litres inspired or expired per minute. However, we can't 
wait for a whole minute, or even several seconds, to calculate it, because we wish to be able to prevent apneas or 
hypopneas lasting even a few seconds, and a PI controller based on an average ventilation over a few seconds would 
be either sluggish or unstable. 

20 [0108] The quantity actually servo-controlled is half the absolute value of the instantaneous respiratory airflow. A 
simple clipped integral controller with no damping works very satisfactorily. The controller gain and maximum output 
ramp up over the first few seconds after putting the mask on. 

[01 09] If we have had a sudden increase in mouth leak, airflow will be nonzero for a long time. A side effect is that 
the ventilation will be falsely measured as well above target, and the amount of servo assistance will be falsely reduced 
25 to zero. To prevent this, to the extent that the fuzzy recent peak jamming index is large, we hold the degree of servo 
assistance at its recent average value, prior to the jamming. 
[01 10] The algorithm for calculating servo swing is as follows: 

error = target ventilation - abs(respiratory airflow) /2 

30 

servo swing = S error * servo gain * sample interval 

clip servo swing to range 0 to 20 cmH 2 0 * iead-in 

35 set recent servo swing = servo swing low pass filtered with a time constant of 25 sec. 

clip servo swing to be at most Jp&w * recent servo swing 

[0111] The instantaneous servo assistance is calculated by multiplying servo swing by the previously calculated 
40 pressure waveform template: 

instantaneous servo assistance = servo swing * n(4>) 
Estimating Instantaneous elastic pressure 

45 

[0112] The instantaneous pressure required to unload the elastic work of inspiring against the user-specified 
elastance is the specified elastance times the instantaneous inspired volume. Unfortunately, calculating instantaneous 
inspired volume simply by integrating respiratory airflow with respect to time does not work in practice for three reasons: 
firstly leaks cause explosive run-away of the integration. Secondly, the integrator is reset at the start of each inspiration, 
so and this point is difficult to detect reliably. Thirdly, and crucially, if the patient is making no efforts, nothing will happen. 
[01 13] Therefore, four separate estimates are made, and a weighted average taken. 

Estimate 1: Exact instantaneous elastic recoil calculated from instantaneous tidal volume, with a correction for sudden 
change in leak 

55 

[01 14] The first estimate is the instantaneous elastic recoil of a specified elastance at the estimated instantaneous 
inspired volume, calculated by multiplying the specified elastance by the integral of a weighted respiratory airflow with 
respect to time, reset to zero if the respiratory phase is expiratory. The respiratory airflow is weighted by the fuzzy rvega- 
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tion of the recent peak jamming index J PEAK , to partly ameliorate an explosive run-away of the Integral during brief peri- 
ods of sudden increase in leak, before the leak detector has had time to adapt to the changing leak. In the case where 
the leak is very steady, J PEAK will be zero, the weighting will be unity, and the Inspired volume will be calculated nor- 
mally and correctly. In the case where the leak increases suddenly, J PEAK will rapidly Increase, the weighting will 
decrease, and although typically the calculated inspired volume will be incorrect, the over-estimation of inspired volume 
will be ameliorated. Calculations are as follows: 

Instantaneous volume = integral of respiratory airflow * (1-J PEAK ) dt 
if phase is expiratory (0.5 < $ < 1.0 revolutions) reset integral to zero 
estimate 1 = instantaneous volume * elastance 
Estimate 2: based on assumption that the tidal volume equals the target tidal volume 

[01 15] The quantity standard swing is the additional pressure modulation amplitude that would unload the specified 
elastance for a breath of a preset target tidal volume. 

target tidal volume = target ventilation / target frequency 

standard swing =_elastance * target tidal volume 

estimate 2 = standard swing * Yl($) 

Estimate 3 ; based on assumption that the tidal volume equals the target tidal volume divided by the observed mean 
respiratory rate RR calcutated previously. 

[0116] 

Estimate 3 = elastance * target ventilation / RR * n(<|>) 
Estimate 4: based on assumption that this breath is much like recent breaths 

[01 17] The instantaneous assistance based on the assumption that the elastic work for this breath is similar to that 
for recent breaths is calculated as follows: 

LP elastic assistance = instantaneous elastic assistance 

low pass filtered with a time constant of 2STD T TOT 

estimate 4 - LP elastic assistance * n(<»/ P BAR 

[01 18] The above algorithm works correctly even if n(<J>) is dynamically changed on-the-fly by the user, from square 
to a smooth or vice versa. For example, if an 8 cmH20 square wave (110^=1) adequately assists the patient, then a 
sawtooth wave (1^^=0.5) will require 16 cmH 2 0 swing to produce the same average assistance. 

Best Estimate of Instantaneous Elastic Recoil Pressure 

[0119] Next, calculate the pressure required to unload a best estimate of the actual elastic recoil pressure based 
on a weighted average of the above. If n($) is set to the smoothest setting, the estimate is based equally on all the 
above estimates of instantaneous elastic recoil. If n(4>) is a square wave, the estimate is based on all the above esti- 
mates except for estimate 1 , because a square wave is maximal at <(>=0 t whereas estimate 1 is zero at <|>=0. Intermediate 
waveforms are handled intermediately. Quantity smoothness runs from zero for a square wave to 1 for a waveform time 
constant of 0.3 or above. 

smoothness = waveform time constant / 0.3 
instantaneous recoil = (smoothness * estimate 1 + estimate 2 + estimate 3 + estimate 4) / (smoothness + 3) 
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[0120] Now add the estimates based on minimum and servo swing, truncate so as not to exceed a maximum swing 
set by the user. Reduce (lead in gradually) if the mask has only just been put on. 

/ = instantaneous minimum assistance + instantaneous servo assistance + instantaneous recoil 

5 

[0121] Truncate I to be less than preset maximum permissible swing instantaneous elastic assistance = / * lead-in 
[0122] This completes the calculation of Instantaneous elastic assistance. 

Desired pressure at sensor 

10 

[0123] 

desired sensor pressure = epap + hose pressure loss + resistive unloading + instantaneous elastic assistance 

15 Servo control of motor speed 

[0124] In the final step, the measured pressure at the sensor is servo-controlled to equal the desired sensor pres- 
sure, using for example a clipped pseudodifferential controller to adjust the motor current. Reference can be made 
toFig. 1 in this regard. 

20 

Device Performance 

[0125] Figs. 21-27 each show an actual 60 second recording displaying an aspect of the second embodiment. All 
recordings are from a normal subject trained to perform the required manoeuvres. Calculated respiratory airflow, mask 
25 pressure, and respiratory phase are calculated using the algorithms disclosed above, output via a serial port, and plot- 
ted digitally. 

[0126] In Figs. 21-26 respiratory airflow is shown as the darker tracing, the vertical scale for flow being ± Usee, 
inspiration upwards. The vertical scale for the pressure (light trace) is 0.2 cmH 2 0. 

[01 27] Fig. 21 is recorded with the servo gain set to 0.1 cmH 2 0/L/sec/sec t which is suitable for subjects with normal 

30 chemoflexes. The subject is breathing well above the minimum ventilation, and a particularly deep breath (sigh) is taken 
at point (a). As is usual, respiratory effort ceases following the sigh, at point (c). The device correctly permits a short 
central apnea (b), as indicated by the device remaining at the end expiratory pressure during the period marked (b). 
Conversely Fig. 22 shows that if there is no preceding deep breath, when efforts cease at (a), the pressure correctly 
continues to cycle, thus preventing any hypoxia. Fig. 23 is recorded with servo gain set high, as would be appropriate 

35 for a subject with abnormally high chemoreflexes such as is typically the case with Cheyne-Stokes breathing. Now 
when effort ceases at arrow (a), pressure continues to cycle and a central apnea is no longer permitted, despite pre- 
ceding deep breathing. This is advantageous for preventing the next cycle of Cheyne-Stokes breathing. 
[0128] The above correct behaviour is also exhibited by a time mode device, but is very different to that of a spon- 
taneous mode bilevel device, or equally of proportional assist ventilation, both of which would foil to cycle after all cen- 

40 tral apneas, regardless of appropriateness. 

[0129] Fig. 24 shows automatically increasing end-inspiratory pressure as the subject makes voluntarily deeper 
inspiratory efforts. The desirable behaviour is in common with PAV, but is different to that of a simple bilevel device, 
which would maintain a constant level of support despite an increased patient requirement, or to a volume cycled 
device, which would actually decrease support at a time of increasing need. 

45 [0130] Fig. 25 is recorded with a somewhat more square waveform selected. This figure shows automatically 
increasing pressure support when the subject voluntarily attempts to resist by stiffening the chest wall at point (a). This 
desirable behaviour is common with PAV and volume cycled devices, with the expectation that PAV cannot selectively 
deliver a squarer waveform. It is distinct from a simple bilevel device which would not augment the level of support with 
increasing need. 

so [0131] Fig. 26 shows that with sudden onset of a severe 1 .4 Usee leak at (a), the flow signal returns to baseline (b) 
within the span of a single breath, and pressure continues to cycle correctly throughout. Although timed mode devices 
can also continue to cycle correctly in the face of sudden changing leak, the are unable to follow the subject's respira- 
tory rate when required (as shown in Fig. 27). Other known bilevel devices and PAV mis-trigger for longer or shorter 
periods following onset of a sudden sever leak, and PAV can deliver greatly excessive pressures under these condi- 

55 tions. 

[0132] Fig. 27 shows an actual 60 second tracing showing respiratory airflow (heavy trace ±1 Usee full scale) and 
respiratory phase as a continuous variable (light trace, 0 to 1 revolution), with high respiratory rate in the left half of the 
trace and low respiratory rate in the right half of the trace. This trace demonstrates that the invention can determine 
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phase as a continuous variable. 

Advantageous aspects of embodiments of the Invention. 

5 Use of phase as a continuous variable. 

[01 33] In the prior art, phase is taken as a categorical variable, with two values: inspiration and expiration. Errors in 
the detection of start of inspiration and start of expiration produce categorical errors in delivered pressure. Conversely, 
here, phase is treated as a continuous variable having values between zero and unity. Thus categorical errors in meas- 
10 urement of phase are avoided. 

Adjustable filter frequency and allowance for phase delay 

[0134] By using a short time constant when the subject is breathing rapidly, and a long time constant when the sub- 
15 ject is breathing slowly, the filter introduces a fixed phase delay which is always a small fraction of a respiratory cycle. 
Thus unnecessary phase delays can be avoided, but cardiogenic artifact can be rejected in subjects who are breathing 
slowly. Furthermore, because phase is treated as a continuous variable, it is possible to largely compensate tor the 
delay in the low pass filter. 

20 Within-breath pressure regulation as a continuous function of respiratory phase. 

[01 35] With all prior art there is an intrusive discontinuous change in pressure, either at the start of inspiration or at 
the start of expiration. Here, the pressure change is continuous, and therefore more comfortable. 
[0136] With proportional assist ventilation, the instantaneous pressure is a function of instantaneous volume into 
25 the breath. This means that a sudden large leak can cause explosive pressure run-away. Here, where instantaneous 
pressure is a function of instantaneous phase rather than tidal volume, this is avoided. 

Between-breath pressure-regulation as a function of average inspiratory duration. 

30 [0137] Average inspiratory duration is easier to calculate in the presence of leak than is tidal volume. By taking 
advantage of a correlation between average inspiratory duration and average tidal volume, it is possible to adjust the 
amplitude of modulation to suit the average tidal volume. 

Provision of a pressure component for unloading turbulent upper airway resistance, and avoiding cardiogenic pressure 
35 instabilities. 

[01 38] Although Younes describes the use of a component of pressure proportional to the square of respiratory air- 
flow to unload the resistance of external apparatus, the resistance of the external apparatus in embodiments of the 
present invention is typically negligible. Conversely, embodiments of the present invention describes two uses for such 
40 a component proportional to the square of respiratory airflow that were not anticipated by Younes. Firstly, sleeping sub- 
jects, and subjects with a blocked nose, have a large resistance proportional to the square of airflow, and a pressure 
component proportional to the square of airflow can be used to unload the anatomical upper airway resistance. Sec- 
ondly, small nonrespiratory airflow components due to heartbeat or other artifact, when squared, produces negligible 
pressure modulation, so that the use of such a component yields relative immunity to such nonrespiratory airflow. 

45 

Smooth transition between spontaneous and controlled breathing 

[01 39] There is a smooth, seamless gradation from flexibly tracking the subject's respiratory pattern during sponta- 
neous breathing well above the target ventilation, to fully controlling the duration, depth, and phase of breathing if the 
so subject is making no efforts, via a transitional period in which the subject can make progressively smaller changes to 
the timing and depth of breathing. A smooth transition avoids categorization errors when ventilation is near but not at 
the desired threshold. The advantage is that the transition from spontaneous to controlled ventilation occurs unobtru- 
sively to the subject. This can be especially important in a subject attempting to go to sleep. A similar smooth transition 
can occur in the reverse direction, as a subject awakens and resumes spontaneous respiratory efforts. 

55 

Claims 

1 . A method for calculating the instantaneous inspired volume of a subject during ventilatory assistance, comprising 
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estimating respiratory airflow, 

calculating an index equal to the fuzzy extent to which any teak Is suddenly changing, and 

calculating instantaneous inspired volume as a function of both the estimated respiratory airflow and said 

index. 

A method as in claim 1 , in which said function is the time integral of the product of said index and said estimated 
respiratory airflow. 
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Fig. la 
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